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Abstract: We build on previous efforts to model CCD sensors, during illumination and collec¬ 
tion of conversions. We use a finite summation of simple, electrostatic field models. The upgraded 
functionality of our framework provides specific predictions for perturbations in pixel boundary 
enclosures (e.g., at the backside window) and the bookkeeping capability to stack those perturba¬ 
tions so that they may be utilized as Greens functions - portable calculation results that may be 
generically applied to a range of precision astronomy related problems that naturally including as¬ 
trometric, photometric and shape transfer issues. We approach the topic of using ancillary pixel 
data, derived from the Greens function and the registered image, to analyze sky data and constrain 
object parameters of astronomical targets. 
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1. Introduction 

A fundamental assumption at work when we apply data from a photo-electronic recording device 
{e.g., CCD) to precision astronomy needs is the notion that the geometric domain tied to each pixel 
can be quantified fo some degree of precision. In fwo previous papers we describe drifl calculafions 
using finite sums of physically mofivated elecfroslafic terms thaf obey boundary conditions OH] and 
subsequent quantitative constraints of the electrostatic parameters using available characterization 
data[0]. 

This paper provides an update on the electrostatic drift calculations using infinite grids and fi¬ 
nite arrays of continuous- and truncated multipoles in arrangements that satisfy electrostatic bound¬ 
ary conditions. In §@ we briefly summarize work already described elsewhere, then in §0 we de¬ 
scribe improvements in modeling and in the modularity of drift calculation results. We ultimately 
(§0) provide some examples for how this framework can be utilized to routinely produce better 
interpretation of laboratory characterization data and how it can in turn be used to significantly 
reduce the systematics floor prior to discovery of any remaining effects, which are typically found 
only after large quantities of high quality exposures are available, acquired only under the best 
conditions. 
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Figure 1: Descriptions of elements contributing to the overall drift field (z) refer to items la¬ 
beled in this figure. Disfribufed dipole momenfs wifh franslafional invariance along symmefry axes 
s are represenfed by fhe exfruded arrows ^cs and ^ck which repeal every pixel; dipole momenfs 
induced by conversions collecled af fhe channel are represenfed by fhe small arrows poinfing down¬ 
ward pij. The backdrop one-dimensional field Ebd is represenfed by fhe arrows poinfing upward 
toward fhe backside eleclrode, which would be 10 pixel unifs away toward fhe top of fhe page. See 
fexl for more delails. 

2. Summary of prior work 

Figure |I] is a perspective represenlalion of a 2x2 pixel section near fhe gafe sfrucfure showing 
fhe elecfroslafic model elemenfs used in fhe drifl calculafions. The disfribufed dipoles fhaf acf as 
baniers againsf conversion loss to adjacenf pixels, are represenfed as exfruded yellow arrows: The 
channel stop dipoles ^cs exfend along fhe parallel address axis [j] while fhe clock dipoles ^ck lie 
in fhe plane of fhe fronl side gale sfrucfure and exfend along fhe serial address axis [i]. When de- 
plefed, fhe p-i- channel stop implanfs form localized, linear bound charge density disfribufions some 
disfance away from fhe clocks: equipofenfial, polysilicon sfrips fhaf pass underneafh fhem in fhis 
picfure. Since fhere can only be a normal elecfric field componenf E = {E -h) -h at fhe surfaces 
of fhose equipofenfials, we model fhe elecfroslafic influence of fhe channel slops using an image 
charge disfribulion to form dipole momenfs ^cs- Similar argumenls are used fo freal fhe influence 
of fhe periodic, allernafing clock polarilies as fhe in-plane dipole momenfs ^ck- The zerolh order 
“backdrop” drifl field Ebd{z) shown is responsible for drawing conversions from approximafely 
100/xm (10 pixels) above fhe plane shown, down fo fhe level of fhe buried channel (OH, §3]). The 
combined effecls from Ebd and fhe fields due to distributed dipole moments ^cs and ^ck - give rise 
to the potential wells in 3 dimensions as well as to saddle point loci encountered by conversions 
that drift along Ebd close to the pixel boundaries (e.g., 0, Fig.l]). Finally, conversions collected 
and confined fo fhe pofenfial wells shown (as franslucenf ellipsoids) - in furn sef up image charge 
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distributions across the polysilicon and an electrostatic influence characterized by pixel specific, 
dipole moments pij. The following list enumerates details of, and deviations from, the nominal 
model described above - in order to match and quantitatively compare with sensor response de¬ 
tails routinely extracted from sensor characterization data. In all cases, pixel boundary distortions 
are computed by successively bisecting adjacent points that connect, via the sum total drift fields 
to neighboring pixels according to methods previously described HH, §4], [0, §2j^]. 

1. The backdrop field. Ebd{z) and its relationship to clock rail voltages and backside bias 
Vbss^ was studied and constrained using available inferential methods. This was described 
in detail earlier [0 §3.1] and comprises a tall pole in the physical description of the device 
under operation. 

2. Edge rolloff. For the edges parallel to the channel stops, the regularly spaced arrangement 
of channel stops is truncated so that the outermost columns are bounded by channel stops, 
but with no more channel stop implants beyond. This abrupt change in the distribution of 
channel stops is sufficient to induce the drop in flat field response with a ~10 pixel scale 
length and ~ 1 pixel astrometric error [0 0 . Presence of a guard ring electrode is modeled 
with a single distributed dipole moment in the plane of the gate structure, and positioned 
directly outside of the last channel stop implant. A negative voltage on the guard ring drain 
bias is modeled as ^gd pointing inward toward the imaging section of the sensor, and results 
in a modest improvement in photo response uniformity, with strongest response improvement 
in a small number of the outermost columns. See Figure 0 top row plots. 

3. Midline redistribution. The bloom stop implant running perpendicular to the channel stops 
is modeled as a single, isolated implant with parameters identical to each of the channel stop 
implants. Its electrostatic influence is assumed to be identical to ^cs hut with its symmetry 
axis parallel to the symmetry axes for ^ck- Variation in the predicted, low light level flat field 
response is induced by adjusting the backside bias Vbss- See Figure 0 bottom row plots. 

4. Dark columns at amplifier boundaries and bamboo. Localized variations in barrier fields 
characterized by ^cSi and/or ^cKj- Isolated [0 Fig.3, tearing onset] and collective [0, Fig.6, 
bamboo] variations in baiiier strengths have been shown to cause pixel geometry distortions 
as influence functions. We suggest that to first order, small lithographic errors that cause 
variations in the separation between parallel clocks result in very little change to the dis¬ 
tributed dipole moments ^cKj, but that errors in the midpoint position, or the location of the 
dipole, result in a first order translation of the pixel boundary between adjacent rows. We 
therefore expect anticorrelation in area errors for pixels within the same column but from 
neighboring rows - this was observed by Smith & Rammer [0] for smoothing rates with 
slope —1 when computing the standard deviation for averaged "lines" (rows). In contrast, 
distributed dipole moments i^cs, should be susceptible to implant errors (ion dosage & depth) 
on top of any lithographic errors. As we showed previously [0 Fig.3], variations in dipole 
strength induces flat field distortion patterns that are different from the pixel boundary jitter, 
and should not necessarily follow the same slope seen in the smoothing rate. 


- 3 - 


5. Brighter-fatter effect. As indicated above and elsewhere [0, 0, these pixel boundary dis¬ 
tortion patterns are induced by the charge pattern already collected midway through the in¬ 
tegration, and influence drift paths of subsequent conversions. The dipole moments pij are 
expected to have magnitudes equal to twice the product of the collected charge in the pixel 
and the depth of the buried channel, measured from the plane of the gate structure. The depth 
of the channel is not known a priori, but can be inferred from details of the photon trans¬ 
fer (mean-variance) curve using a set of predictive calculations. Our previous, strawman 
attempt to do this using a simplistic, nearest-neighbor only boundary distortion model OJ, 
§4.5, Fig.7] would be superseded with our current model, because detailed pixel area dis¬ 
tortions we model now carry enough information to adhere to basic conservation rules. See 
Figures ^ (upper left panel) and ^ (top row panels) for examples of how pixel boundary shifts 
are calculated in response to various levels of charge integration levels and its distribution. 

6. PRNU. In the course of computing pixel geometry response (at the backside surface) due to 
small geometric and electrostatic deviations from the regular, repeated grid depicted in Fig¬ 
ure |T], we realized above (item #0) that the influence functions on pixel geometry due to litho¬ 
graphic or process geometry differences are significantly different in scale and symmetry, 
and these could potentially explain the marked difference in photon response non-uniformity 
(PRNU) smoothing rates 0] and power spectra density slopes [0] along the parallel and serial 
address axes. We haven’t any more details on this subject now, but plan to investigate further 
in the near future, because differences in expected area structure functions along the two 
address axes could be useful in removing degeneracy between pixel area and pixel quantum 
efficiency in any detailed decomposition of sensor flat field response. 

7. Charge transfer efficiency. None of the work described here pertains to CTE, but we antic¬ 
ipate needs in the future to properly simulate the effects of traps unexposed by the fat zero 
provided by sky background, that happen to have de-trapping times comparable to parallel 
transfer timescales. There is already evidence |0] that would suggest localized fluctuations in 
parallel transfer efficiency at medium light levels, and is perhaps not reflected in the extended 
pixel edge response (EPER) method of quantifying the sensor’s CTE for device acceptance 
purposes. Such CTE mechanisms may have non-negligible impacts on PSE estimation for 
high signal-to-noise images and have yet to be quantified. 

3. Improvements to the current modeling framework and recalibration 

In §0 we briefly described the geometry where we apply the individual electrostatic contributions to 
the total drift field {cf. OH, §§3-4]). Since publication of the second paper |0], several important 
changes have been implemented, which are listed here: 

1. Dipole moment ratio calibration against observed flat field covariance ratios. Relative 
values of distributed dipole moments ^cs and ^ck that nominally reproduce observed ratios 
in covariance terms AlnAoi/AlnAio ~ 3 [0, Eig.7], must be recomputed because of the new, 
in-plane arrangement of dipole moments ^ck with the alternating polarity shown in Eigure 
Because of the alternating polarity shown, it could be argued that the barrier influence of 
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the nearest two clock pairs may resemble that of a quadrupole moment instead of a dipole 
moment - so a different radial dependence would be in effect in the far field approximation. 
The updated parameter determination is shown in Figure 0 and discussed at more depth in 
§0- 


2. Update to the method of images strategy. Another side effect of the new dipole mo¬ 
ment orientations for is that the plane of the gate structure is no longer be treated as an 
equipotential plane imparting its boundary condition on the electrostatic solution |]I], Eqs.4.10 
& 4.12]. The modification of those equations that is consistent with a single equipotential 
plane at the backside entrance window (for field confribufors fhaf depend on ^ck) is to in¬ 
clude only a single image for each charge configuration contributing to in place of the 
expansion given in those equations. Equation 4.12 would be altered to truncate the sum from 
0 < £ < 1. An updated representation of Equation 4.10 now includes a mapping of the image 
dipole moment: 

_ i ^ 2[ff • 
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3. Influence of an inflnite, periodic arrangement of barrier terms. The regularly spaced, 
electrostatic model contributors, summed over the infinite grid - is now evaluated using the 
following periodic expression we derived analytically: 
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where s is the symmetry axis of the distributed dipole and k is the axis perpendicular to the 
plane of the grid of distributed dipoles, which are in turn spaced at regular intervals with 
(pixel) spacing p along the direction ei_ in the gate structure plane. It follows that (k,5,ej^) 
form an orthonormal set of unit vectors. Use of this closed expression, together with a single 
image term per electrostatic elements ^ck rather than a finite expansion, greatly accelerates 
evaluation of the drift field wifhin the pixel volume. Any deviation from the perfectly uniform 
and regular arrangement in distributed dipoles is easily be implemented by cancelation and 
modification of individual dipole terms in the expression, since the principle of superposition 
applies. 
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4. Discontinuities in the distributed dipoles. Expressions were also derived to represent the 
influence of distributed dipole moments that are truncated - that feature finite edge effect, 
field components parallel to the symmetry axis s nearby the truncation coordinate. These 
components are useful in modeling the effects of tearing - where apparently the hole con¬ 
centration within the channel stops can abruptly transition from zero {i.e., depleted) to some 
saturation level where transfer efficiencies drop nearly to zero. The superposition of such 
truncated, distributed dipole moments in adjacent columns can reproduce the sort of bimodal 
features seen in flat field response near the position of the change in hole concentration: this 
is the contour seen in flat field data where “tearing” is a problem {cf. lU, Fig.5]). It also 
applies to the “bamboo” phenomenon {cf. |0, §4.2]) for reasons described there. 

The expressions derived for the truncated, distributed dipole moment are: 

Aii£Qesir\5E (x|xo, f) = (2 — 3 cos do + cos^ 0o) [I • ^_l] r± 

+ sin^ do ■r±]s- (1 - cos 60)!, (3.10) 

cos00 ^ s-r. (3.11) 

Here, the distributed axis of symmetry unit vector s points in the direction of the truncation, 
i.e. from the truncation point xo away from where the distributed charges are situated. As 
required, this expression converges to the familiar distributed dipole expression in the lim¬ 
iting cases. However, the truncated distributed dipole model does not lend well to forming 
concise expressions such as the infinite sum above, so a summation over a finite number of 
electrostatic components is performed when necessary for certain fixed pattern features {e.g., 
discontinuous features associated with tearing onset []I], Fig.3, bottom row plots]; “bamboo” 
features 0, Fig.6]). 

5. Porting drift calculation results to modeling platforms. We have recently devised a way 
to make drift calculation results portable so that contributions from various effects can be ef¬ 
ficiently scaled and combined in downstream simulation programs - where pixel boundaries 
are treated as simple polygons with N ordered vertices {xk,yf) with k = d..N — \, closing 
them with {x^^Jn) = (-^oAo)- In the limit that drift calculations produce null results, the 
pixel boundaries then follow the usual square contours, consistent with the fundamental as¬ 
sumption of common pixel data pipeline routines. The information caiTied from the drift 
calculations include, for each segment of the pixel boundary, specific to a given aggressor 
offset vector and amplitude: 

• Fateral, two dimensional displacements of the positions at the backside window that 
map to the corners of the registered pixel for that segment (4 numbers); 

• An array of normal displacements perpendicular to the segment. The number of sam¬ 
ples, together with an array of the displacement values with uniform spacing along the 
nominal segment axis. 

Our downstream simulator loads these distortion patterns as templates, scales and superposes 
the effects in according to rules of the model, and efficiently evaluates the pixel boundaries. 
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Boundaries represented as ordered polygon vertices are made accessible from the data struc¬ 
tures or classes that represent each pixel. The polygonal pixel structure (together with higher 
level summary quantities) lends to efficient point-in-polygon algorithms to test domain mem¬ 
bership, and to robust, single pass, geometric formulae. The following pixel moment expres¬ 
sions were derived by generalizing the shoelace formula for polygon area: 
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The sign of Eq. [3.12| corresponds to a specific choice of chiralify for the polygonal vertex 
list. The quantities above are used to evaluate distortions to pixel area (5lnA,y), pixel astro- 
metric shift vectors {e.g., pij -x = [Ix,yA,y]/A,y) and pixel ellipticities {5eijj = [IxXijAij — 
lYYijAij]/[lxXijAij-ylYYijAij]-, de2,ij = 2IXY,yA,y/[Ixx,yA,y-hlYY,yA,y]). Examples for 
calculations of pixels in the template are shown in Eigure 0 and those evaluated at the end 
of an exposure for a bright, stellar source for three specific operating conditions are given in 
Eigure 

This internally consistent partition model may have wide applicability to detailed modeling, 
starting with validation against available laboratory characterization data: flux dependence of flat 
fields and the shift-covariance statistics they express, modulation transfer function (MTE) measure¬ 
ments, and spot projector illumination using well characterized beams. Proper validation should 
include quantitative comparisons of the flux dependences for the above tests. In our opinion, the 
most fundamental property and largest signal to reproduce is the shape of the mean-variance re¬ 
lation obtained from photon transfer measurements. Once validated, the same modeling platform 
may be used to compare and quantify differences between incident flux distributions and recorded 
conversion distributions. 

4. Examples 

The remainder of this paper is devoted to depicting some results made possible by the updated 
modeling capabilities described above. In the future, we plan to identify better verification metrics 
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against sensor characterization data acquired under a fixed set of operating conditions - low light 
levels and voltage settings. Spot projector data [@] is also of extreme interest, even if the incident 
beam is not well understood: we are equipped to model and compare changes in the distribution of 
recorded images as a function of integration level. 

4.1 Tuning dipole strength ratio for consistency with observed covariance term ratio 

To tune the dipole strength ratio, we started with the assumption that the bloom stop implant ^bs 
is identical to each of the channel stop implants (modulo its symmetry axis s orientation). We 
matched the low-flux flat field response signature to the computed “fixed pattern” pixel areas there. 
Using a backdrop field profile Ebd{z) corresponding to the Vbss = —50V setting for that data, 
we converged on placeholder for the implant distributed dipole moment = —14.3i^o^\ where 
as previously, i§o = lO^^qe = 1.6 x lO^^^Coul. This choice also provided edge rolloff pixel area 
profiles that matched reasonably well when compared against the flat field response signature there. 
Considering the value of ^cs to be fixed, the value of the clock rail induced dipole moment ^ck 
was constrained by computing the ratio of geometric area variations for neighboring pixels due to 
presence of accumulated charge in the “aggressor” pixel^. This ratio is plotted in Figure ^ (upper 
left panel) for a range of backside bias voltages and for a selection of trial clock dipole values. 
We identified a “fiducial performance” region of this parameter space that exhibits the observed 
ratios in geometric area variations, inferred directly from covariance ratios [01. The “tuned” value 
for the clock dipole magnitude was then chosen to be \^ck\ = 5i^o- This specific combination of 
contributors to now features pixel neighbor area perturbations that respond to collected charge 
(the brighter-fatter effect) with proper proportions. It should be noted that a joint fit could also 
have been performed to match the various observables. For example, such fits may favor a different 
dipole strength ^bs for the bloom stop bander, a different acceptor density Na for the photosensitive 
bulk, or a different field strength offset when Vbss = 0 and the device is fully depleted. We have 
not yet performed such fits that may further tune parameters of this model. 

The influence of collected charge at the channel remains to be calibrated. This is a straight¬ 
forward step that is best done using constraints of the mean-variance relation. It is equivalent to 
determining the depth of the buried channel in this simplistic model, because the dipole moments 
Pij scale with that dimension. We also expect that the depth or position of the potential well will 
vary with the number of conversions collected. These details will depend on a number of factors, 
including the doping profile of the n-type channel and the attractive influence of the image charge 
that, combined together with the collected charge cloud in the channel and the presence of ad¬ 
jacent, similar, conversion-induced dipole moments in neighboring pixels, self-consistently form 
the dipole moments pij. We assume that such details may be most easily accessed from photon 

'This dipole strength is approximately twice as large as that determined previously [0, Fig.5]. This difference is 
partially explained by the change in the symmetry of the electrostatic problem, where previously two equipotential 
planes were used, and a truncated infinite expansion of image charges were used. 

^For convenience and for an initial tuning this model, we adjust the value of ^k- What is really being adjusted in so 
doing is the position within the photosensitive bulk where the backdrop field Ebd(z) is nearly canceled by the collective 
5E^. There are multiple solutions that could provide this knob, including a modification of the assumed doping profile 
through the photosensitive bulk that provide the built-in field [HI, 0]. We have not identified other independent observables 
that constrain the built-in field, so for now we choose to leave Ebd{z) fixed. In the future, the dependence of covariance 
terms on BSS may be used to jointly constrain Ebd{z) and SE^ . 



transfer data sets. For now, we consider only potential well depths that are constant with integrated 
conversions. 

The resulting brighter-fatter template was generated for pij = — 1.05,o5,/po^ (where po = 
lO^qeftm) and is available in the form of pixel specific boundary contours (Figure ||, upper left 
panel) and the pixel geometric moments obtained using the shoelace formulae given above. The 
latter are plotted in Figure upper right and lower panels. The overall radial dependence of these 
parameters (lower panel) provides an overview of the brighter-fatter effect, with radial dependences 
evaluated by fitting over the range 4 < (A/)^ -|- (Ay)^ < 25pix^. The following radial scaling de¬ 
pendences were obtained for the fitting range: A5lnA,y/Alnr^ « —1.288, A5ln(l —Pr,ij) /Alnr^ Ri 

— 1.033 and A5ln(l — er,y)/Alnr^ « —1.410. Here, the radial component of ellipticity for pixel 
{ij) is defined as Erjj = {hr — he)/{Ixx + lyy)- Notice that these slopes are dependent on the fit 
domain applied. Two features should be pointed out: First, the geometric distortions for distances 
within the fit domain and farther - should be relatively insensitive to the influence of ^cs and ^ck, 
because most of the boundary shifts occur as a result of curved drift lines where the backdrop field 
Ebd is smallest. This also implies that this component of the geometric distortions should have 
some chromatic dependence. The geometric distortions appear to fall off rapidly beyond 10 pixels 

- which happens to be equal to the assumed thickness of the sensor. This is reassuring, because 
influences fall off more quickly than and a limited number corrections might be required in 
any pixel pipeline that utilizes them. Second, the large deviations from the power law fits for the 
nearest pixels are unique and the short distance effect is clearly not isotropic (a consequence of the 
tuned dipole moment ratios). The drift contribution to these geometric distortions must be small, 
because the lateral electrostatic influence of a dipole field for points along the dipole’s axis van¬ 
ishes. What appears to govern the pixel boundaries closest to the collected conversions, are the 
distortions of the saddle point loci experienced by charges flowing along the drift lines close to the 
pixel boundaries. This is the component of pixel geometric distortion that can’t be simulated with 
simplified elecfrostatic drift models. Because these shifts in pixel boundaries occur close to the 
gates and channel stops, this contribution should be nearly achromatic. 

4.2 Large amplitude fixed pattern features 

We evaluate the pixel boundaries near the edges (parallel to the parallel address axis) and those 
crossing the midline bloom stop. Using the fiducial performance paramefers outlined above, we 
obtain the pixel contours, pixel areas A,y, pixel position shifts pij and pixel ellipticities eij for the 
affected pixels. They are summarized in Figure where we have overlaid those quantities for 
different guard ring drain bias induced dipole moments (for the edge response case). We have 
also overlaid quantities for different backside bias values (for the midline redistribution case). In 
all cases we predict subtle variations in the detailed pixel shapes when evaluated at the backside 
surface - the edges often appear to be scalloped in shape and the notion of parallel pixel boundaries 
is generally violated. 

4.3 Flat field simulations that include dynamic pixel boundary distortions 

Here we simulate flat fields to predict the relevant statistical properties of those illuminations. Re¬ 
call that in above, the clock dipole strength fcx was adjusted so that the ratio of pixel area 
perturbations would match the measured pixel covariance ratio: 5lnAio/5lnAio = Cio/Coi ~ 3. 
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Figure 2: Edge rolloff response (top plots) and midline redistribution (bottom plots) large amplitude 
fixed pattern feature pixel calculations. For each effect, pixel boundary contours are shown for two 
rows or columns (left hand plots), while geometric parameters A,y, pij and Eij are shown (right hand 
plots). Predictions are given for adjacent guard ring drain bias settings (top) and for backside bias 
settings (bottom). 


This is a check to test the validity of that assumption. Starting with the brighter-fatter template de¬ 
scribed in Figure 0, we populate a stack of 1024 small images with a randomly distributed, constant 
flux to represent a flat field illumination. This is done in steps of on average 1000 conversions per 
pixel, using a Poisson deviate together with the pixels’ current areas Aij. Between population steps, 
some statistical calculations are performed on the recorded charge patterns in the image stack, and 
those patterns are used to reevaluate the pixel areas that feed into the next population cycle. The 
statistical calculations are performed much in the same way as laboratory data would be: for every 
unique pair of images, a difference image is computed. For specific serial and parallel shift (lag) 
values, product of the difference values is summed to contribute to an overall covariance. The 
grand total sum is divided by the total number of electrons contributing to the calculation (2/rNpix) 
to yield the covariance value for lag (A/, Ay). When performed for a single image pair, the 
covariance is simply the numerator in the expression for a difference image’s autocorrelation ele¬ 
ment. 

The calculation is complete after 250 such population and analysis cycles. Two assumed 
responsivities were applied to the mechanism. The first assumes that the brighter-fatter template 
corresponds to an isolated channel population of 10^ conversions, at which point the central pixel 
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area shrinks to Aqo = 0.9033; adjacent pixel areas grow modestly to become Aqi = 1.0077, Aio = 
1.0027 and All = 1-0033 (5lnAoo = —1 x 10^^/e^). The second assumed that the influence is 3.3 
times as strong, for the same channel population (5lnAoo = —3.3 x 10^^/e^). When we compare 
these computed Cqo curves to carefully prepared photon transfer curves ([0, Fig.4]), it is clear that 
a more realistic responsivity of this effect would be about a factor of 3 smaller than the first case, 
or 5lnAoo = —3 x 10^®/e^. 

Figure ^ provides the results of these calculations, where we plot only the mean-variance re¬ 
lation Coo and the three covariance terms Coi, Cio and C\\. One plot is given for each responsivity 
assumed. Notice that the transfer of variance from the central pixel containing statistical fluctu¬ 
ations, out to adjoining pixels whose areas adjust in response, is complex enough that pixel area 
distortion ratios may not be accurately inferred directly from covariance term ratios, at least for 
the calculation provided here. Potential complications to the working model include at least the 
following: 1) only a fraction (~ 35%) of the a given pixel’s area perturbation due to its contained 
charge excess or deficit, is recovered by the 8 that directly adjoin it; 2) that there are multiple 
aggressors that affect an individual pixel’s area perturbation; and 3) this generalized photon trans¬ 
fer calculation exhibits linear behavior only for the lowest light levels, where these calculations 
clearly fall short statistically. For the near term, we would like to quantitatively understand details 
of the mean-variance curve’s shape (Cqo) in sensor characterization data, since to first order its area 
perturbation is largely driven by its own charge excess or deficit. 

4.4 Point source response using the brighter-fatter template as a Greens function 

An initial investigation into how drift calculations can be used to predict errors in PSF estimation, 
we utilized the same brighter-fatter template in a Monte-Carlo to simulate image formation for 
an unresolved point source. For the PSF model, we drew sample positions one at a time from an 
ellipsoidal Gaussian function, with FWHM set by two possible seeing values (0.6" and 0.8") and 
an ellipticity of 0.05. The major axis is rotated 30° counter clockwise from the serial address axis, 
with centroid positioned at (0.25,0.25) pixel units measured relative to the central pixel’s centroid. 
Starting with an undistorted array of pixel boundaries, sample positions were partitioned into the 
grid of registering pixels using a point-in-polygon algorithm. Conversions were accumulated 1000 
at a time, which typically boosted the central pixel content by about either 80 or 50 conversions, 
depending on the seeing level in effect. Moments of the accumulating image were computed, as 
were changes in the image since last update. Following this, the brighter-fatter template was used to 
update the pixel boundaries based on the current, collected charge pattern. This operation closely 
resembles what was done for the flat field simulations {cf. §|4.3D, except that conversions were 
added one at a time. 

Results of this image accumulation calculation are given in Figure ^ As in the case of the flat 
field integration examples above, the same two brighter-fatter responsivities were applied, which 
correspond to central pixel shrinking rates of 5lnAoo = —1.0 x 10^^/e^ and 5lnAoo = —3.3 x 
10^^/e^. A total of three calculations are summarized. For 5lnAoo = —1.0 x 10^^/e^, seeing 
levels of 0.6" and 0.8" were simulated, whereas for 5lnAoo = —3.3 x 10^^/e^, only the 0.6" 
FWHM seeing was simulated. The figure shows the computed central pixel value as a function of 
source counts, measures of broadening arranged both by source counts and central pixel counts, and 
properties of the recorded source ellipticity arranged by central pixel counts. The latter two show 
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the systematic changes in PSF parameters (recorded source ellipticity and recorded source rotation 
angle) that arise from the pixel boundaries that shift in response to the ongoing accumulation of 
conversions. Not surprisingly, the systematics appear to depend on source integration level, seeing 
and magnitude of the brighter-fatter effect. 

If such a brighter-fatter model is validated adequately against available lab data, it is natural 
to expect that the geometric distortions predicted on a pixel-by-pixel basis can be used in pixel 
data pipelines to partially cancel the effects of the pixel boundary shifts that apparently occur while 
the conversions accumulate. The recorded image would be used to determine the pixel boundaries 
at the end of the accumulation using the brighter-fatter template as the Greens function, whereas 
initial pixel boundaries are assumed to be undistorted. The trajectory connecting the two bound¬ 
ary configurations then depend on details of the flux distribution relative to the undistorted pixel 
boundary grid and its time history. Figure ^ shows the final configuration for a 15 second illumina¬ 
tion by a bright star in r band, AB~15 magnitudes for the three conditions described above. The 
maps show the distorted pixel boundary contours, the pixel area distortions pixel astromet¬ 
ric distortion and pixel ellipticity arranged by pixel positions of the recorded image. 
Notice that accurate background subtraction of the image should also depend to some degree on 
^fijiai ^g ^ range of choices as to how these maps may be generated and used in pixel data 

pipelines, and we leave those details for a later discussion. 

The unresolved source simulation described here also provides access to instantaneous pixel 
accumulation rates. Since tracking errors were not included, the only variations in the rates are 
due to pixel geometric distortions mediated by the brighter-fatter template. Figure 0 shows this for 
select pixel rate averages in the neighborhood centered on the star. In some extreme cases (e.g., 
0.6" FWHM, 5lnAoo = -3.3 x 10^PSF counts 2 X lO^e ) rates at outer pixels overtake 
the central pixel’s rate. In principle, small tracking errors can also lead to image elongations and 
astrometric biases that scale with source brightness - and these errors should be different along the 
two pixel address coordinates. Studies that address these effects can be performed with a modest 
extension of the existing framework. 


5. Conclusions 

We have introduced a quantitative modeling framework that can be used to affect a widespread 
departure from the approximately correct assumptions that pixel boundaries remain, over the dy¬ 
namic ranges of interest, fixed, independent of sensor position and local illumination history. Rich 
datasets can be used to tune parameters of this drift model and be tested against validation re¬ 
quirements. Results of the tuned model are portable and can be applied in superposition outside 
of the sensor drift and partition calculator. We consider how the overall partition model may be 
broken up into influence components that act on longer length scales (and should have moderate 
chromatic and backside bias voltage dependences), and influence components that act on shorter 
length scales, that have significantly reduced dependences on those variables. The superposition of 
brighter-fatter templates on top of effects that inject bias {e.g., edge rollof, tracking errors) - should 
produce brightness-correlated astrometric and shape transfer errors that should be examined. 
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Figure 3: Properties of the brighter-fatter pixel distortion template. Upper left: the dependence of 
the area distortion ratio 5lnAoi/5lnAio on the backside bias BSS for a range of clock induced 
dipole moment values ^ck- The round black dot indicates the fiducial performance (parameter 
choice) for the remainder of this paper. Upper right: a false color representation of the pixel 
geometric parameters in the template, arranged by address: areas A,/ (top two), pixel astrometric 
shifts pij (middle two), and pixel ellipticity Sij (bottom two). Bottom: pixel geometric parameters, 
this time plotted against separation (r^ = f + j^) from the central pixel. Of these parameters, only 
5lnA,y are constrained directly by high signal-to-noise, exposure corrected flat field images. The 
radial dependence agrees wifhin a few percenf of fhe dependence found for a complefely differenf 
device (DECam) where a full season’s worfh of r-band dome flafs were joinfly used for 59 devices 
(ins Fig.4]) - buf fhis driff model also explains fhe nafure of fhe high significance, high amplifude 
oufliers seen fhere. The assumed dipole of fhe aggressor for fhis femplafe was poo = —Po^- See 
fexf for more information. 
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Figure 4: Covariances Cij of difference images where a lag (/, j) has been applied in the calculation. 
The covariance Cqo is simply the variance. Other terms Cqi, Ciq and C\\ are numerators in the 
autocorrelation of the difference images. These terms rise as for low light levels, where also the 
variance deviates only linearly from the Poisson behavior linearly with jx. At intermediate values 
of jx, where the deviations from Poisson behavior become dominant, Cij grow linearly with jx, and 
then flatten off at larger values. The straight lines showed here provide guides for Poisson behavior 
in Coo, and for the linear growth region for Coi and Cio, based on area perturbations 5lnAoi and 
5 In A 10 computed from the brighter-fatter template. Two different response levels for the brighter- 
fatter effect were assumed: 5lnAoo = —1 x 10^^/e^ (left) and —3.3 x 10^^/e^(right). Carefully 
measured curvature of the mean-variance curve for specific operating conditions [0, Fig.4] would 
calibrate the response to ~ —3 x 10^®/e^. Pixel areas were reevaluated for every 1000 electron 
increment to the flat field level, or once per plotting symbol shown. 
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Figure 5: Summary of the three focused image simulations performed for pixel recording prop¬ 
erties. Upper left: a display of the pixel boundary contours for the brighter-fatter template used. 
Upper right: central pixel counts, arranged by PSF counts. Middle left: recorded PSF area a? -|- oj, 
arranged by source counts. Middle right: FWHM contribution by the brighter-fatter effect (RSS 
difference) arranged by central pixel counts. Lower left: ellipticity of the recorded image, arranged 
by central pixel counts. The drop in ellipticity occurs because the input source had a modest ellip¬ 
ticity (0.05) oriented at 30°, whereas the brighter-fatter effect is anisotropic. Lower right: recorded 
source orientation, based on ellipticity components. This shows the anisotropic properties of the 
brighter-fatter effect, as brighter objects will be systematically rotated toward a 90° orientation. 
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Figure 6: Pixel specific geometric parameters, corresponding to the final, recorded image. Left col¬ 
umn: 0.6" seeing, 5lnAoo = —1 x 10^^/e^; center column: O.S" seeing, 5lnAoo = —1 x 10^^/e^; 
left column: 0.6" seeing, 5lnAoo = —3.3 x 10^^/e^. Each column gives the final pixel boundary 
contours (first row), final pixel area (second row), final pixel astrometric shifts (third and fourth 
rows), and final pixel ellipticities (fifth and sixth rows). Corrections to individual pixel values 
may be performed using values from these maps to perform better background subtraction and to 
estimate source parameters in a way that is less subject to systematics of the brighter-fatter effect. 
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Figure 7: Pixel specific conversion arrival rates (for specific pixels distributed about the PSF cen¬ 
troid), arranged by PSF counts. Top: O.b” seeing, 5lnAoo = —1 x 10^^/e^; middle: O.S" seeing, 
5lnAoo = —1 X 10^^/e^; bottom: O.b" seeing, 5lnAoo = —3.3 x 10^^/e^. See text for more 
information. 
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